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Abstract. We show that results of a simple dynamical gedanken experiment inter- 
preted according to standard Newton's gravitational theory, may reveal that three- 
dimensional space is curved. The experiment may be used to reconstruct the curved 
geometry of space, i.e. its non-Euclidean metric ^Qik- The perihelion of Mercury advan- 
ce and the light bending calculated from the Poisson equation ^g'^^ViVk^ = —AnGp 
and the equation of motion = ma^ in the curved geometry '^gik have the correct 
(observed) values. Independently, we also show that Newtonian gravity theory may 
be enhanced to incorporate the curvature of three dimensional space by adding an 
extra equation which links the Ricci scalar with the density of matter p. Like 
in Einstein's general relativity, matter is the source of curvature. In the spherically 
symmetric (vacuum) case, the metric of space '^gik that follows from this extra equation 
agrees, to the expected accuracy, with the metric measured by the Newtonian gedanken 
experiment mentioned above. 
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1. Introduction 

Newton's theory of gravity was formulated in a fiat, Euclidean 3-D space but its basic 
laws, i.e. the Poisson equation and the equation of motion. 



make perfect sense in a 3-D space with an arbitrary geometry ^gik- Indeed, the curvature 
of space is (potentially) present in Newton's theory. It is easy to argue that the 
"centrifugal" acceleration of a particle moving with velocity \^ on a circular orbit 
equals = V"^ /TZ, and the "gravitational" acceleration in the gravitational field of 
a spherically symmetric body with the mass M. equals = GM./f'^, with TZ being 
the curvature radius of the circle, and f being its circumferential radius. In fiat, i.e. 
Euclidean, 3-D space these two radii are equal, TZ = r, but in a space with a non-zero 
gaussian curvature Q, they are different, TZ ^ f. Therefore, by measuring centrifugal 
and gravitational accelerations one may independently measure TZ and f, and thus 
experimentally find whether the space is fiat (Euclidean) or it has a non-zero Gaussian 
curvature ^ 7^ 0. Based on that, Abramowicz has recently suggested in [1] that a 
Newtonian physicist could experimentally determine the metric ^gik of the real physical 
3-D space and calculate, according to ([I]) and ([2]), the perihelion of Mercury advance 
and the light bending effects. In this paper we follow this suggestion and calculate both 
effects within Newton's theory. Surprisingly, the values of the perihelion advance and 
the light bending agree (to the expected order of M/r) with predictions of Einstein's 
theory. Here M is the "geometrical" mass of the spherical gravitating body expressed in 
the convenient "geometrical" units G = 1 = c. It is connected to the mass M. expressed 
in the standard units by M = GAi / (? and has the dimension of length. 

Another point discussed in this paper is based on the following two remarks: 
{%) Obviously, Newton's gravity theory is a limit of Einstein's general theory of relativity. 
Should the limit necessarily correspond to ^ = 0? Perhaps not, because Newtonian 
physicists could discover within Newton's theory that G ^ 0. (M)They could also 
discover that the curvature of space depends on the distance from the gravity center, 
Q = ^(r). This would suggest to them, again within the framework of Newton's theory, 
that gravity and curvature are not independent, but instead they are somehow linked. 
Here we suggest that it is possible to establish the link within an "enhanced" version 
of Newton's theory, by adding to its standard version defined by ([T]) and ([2]) an extra 
equation. 



where is the Ricci scalar corresponding to ^gtk, p is the density of matter, and k 
is a constant. Equations ([1]), ([2]) and ([3]) define our enhanced version of Newtonian 
gravitational theory. In the special case of a spherically symmetric, vacuum (p = 0) 



Fi = mai, 



AnGp, 



(1) 
(2) 



= 2kp, 



(3) 
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Figure 1. For a circle placed in a curved space (here on a curved 2-D surface), its 
geodesic radius r,, circumferential radius f, and curvature radius TZ are all different, 

space, they uniquely lead to the 3-D metric of the form, 

ds' = (l - ^) dr' + {de^ + sin^ ed(j)') , (4) 

where tq is a constant. A choice tq = 4 M leads to correct values for both the perihelion 
advance and the light bending effect^. 



2. The three radii of a circle 

Consider a two dimensional curved, axisymmetric surface with the metric 

ds^ = drl + [f{r,)fd(p\ (5) 

and a family of concentric circles r* = const in it. One of them is shown in Figure [2j 
Obviously, is the geodesic radius and f is the circumferential radius of these circles, 

(geodesic radius) = / ds\^=coryst = (6) 

Jo 

1 f^"" 

(circumferential radius) = — / ds\r =const = (7) 

27r Jo 

Let r* = r~^(50 be a unit vector tangent to the circle. From the Frenet formula, 
dr* 1 

— — = — — A\ where A* = (unit normal to the circle), (8) 

(JjO / V' 

one deduces that the curvature radius TZ may be defined by. 



(curvature radius) 



dr* dTi 
ds ds 



-1/2 



n. (9) 



I Assuming that light moves along geodesic lines in space. 
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Two useful formulae for the curvature of the circle, /C = 1/7^, and for the Gaussian 
curvature Q of the surface with the metric ([5]) read, 

« = (10) 



Formula (fTOj) follows from (|9]). For derivation of (fTT!) see e.g. [2j, Section 3.4. 
3. Equations of motion 

Let us consider a curve in space given by a parametric equation, 

x' = x\s), (12) 

where are coordinates in space, and s is the length along the curve. If a body moves 
along this curve, its velocity equals, 

, dx* dsdx"- 

w = — = ^ — — = vt\ (13) 
dt ds ^ ^ 

Here v = ds/dt is the speed of the body and r* = dx^ /ds is a unit vector tangent to the 

curve (fT2l) . i.e. the direction of motion. The acceleration may be calculated as follow, 

• dv'' dsdivT') 2 f dr'\ ,• dv ^. 
a' = — = - — = v {^] + t'v—. (14) 
dt dt ds \ds / ds 

Assuming circular motion with constant velocity, v = const, and applying ([9]) to 
calculate the term in brackets, we arrive at 

a' = v^—X\ (15) 

which is the well known formula for the centrifugal acceleration. 

Consider now circular motion around a spherically symmetric center of gravity. The 
Newtonian equation of motion, = ma*, takes the form, 

_V*$ = ^2^A\ (16) 

where = — mV*$ is the gravitational force, and $ is the gravitational potential. 
Three quantities characterize motion on a particular circular orbit: the angular velocity 
Q, the angular speed v, and the specific angular momentum C They are related by, 

V = m, (17) 

C = fv = r^n. (18) 



Using ( ITSl) . and multiplying its left side by Aj, we transform the equation of motion (fT6|) 
into a form which will be convenient later, 

A.V*$ = 4^- (19) 

In this expression. A* is a unit, outside pointing, vector. Here "outside" has the absolute 
meaning — outside the center, in the direction towards infinity. We will calculate the 
left-hand side of this equation in the next Section. 
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4. Newton's gravity and Kepler's law 

In an empty space, the gravitational potential $ obeys the Laplace equation, 

Vi(V*<l>) = 0. (20) 

Let us integrate (l20l) over the volume V that is contained between two spheres, 
concentric with the gravity center, with sphere $i being inside sphere $2- We transform 
the volume integral into a surface integral, using the Gauss theorem 

0= / Vi{V'<i>)dv = [ (y'<i>)Ni^^ds+ [ (y'<^)NPds. (21) 

JSI J%2 

The oriented surface elements on the spherical surfaces $i and S2 may be written, 
respectively, as 

Nf>dS = -X.dS, N^^^dS = +XidS, (22) 

therefore, 

f {V'^)XidS= I {y'^)\id^. (23) 

This means that the value of the integral is the same, say So, for all spheres around 
the gravity center. In addition, because of the spherical symmetry of the potential, the 
quantity (V*$)Aj is constant over the sphere of integration. Thus, 

So = (V'$)A, [dS = Attt' (V^<I>)A„ (V*<I>)A, = = (24) 

Combining with f lT9|) . we may finally write, 

£2 = GMn. (25) 

This is the Kepler Third Law. Using natural units for radius and frequency, 

CM r-^ r 

R^^ — ^M, n,- — ^-, (26) 

we may write the formula for the Keplerian angular velocity as, 

5. Epicyclic oscillations, the perihelion advance 

Suppose that we slightly perturb a test-body on a circular orbit. This means that its 
angular momentum will not correspond to the Keplerian one, C^, given by ( 125|) . but will 
be slightly different + 6C^. There will be also a small radial motion with velocity 
(5r*) and acceleration {6r^:). From f|T9|) it follows that 
GM £2 ^ 

= ^2^^ 

Keeping the first order term in equation (1281) . and using 

dC^ 

6C' = ^{6r,), (29) 
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we arrive at the simple harmonic oscillator equation, 

uj\6r,) + {6r,) = 0, (30) 
where uj is the radial epicyclic frequency, 

(31) 



Using equations f l2^ and f l2B]) . we may write the expression for the epicyclic frequency 
in the form, 



or comparing this with flTTI) . 



(32) 



C^V _ ( ^\ _ ( dxV fcPf\ 

In a fiat space, r^, = r = 71, and therefore u = Q, which implies that the slightly non- 
circular orbit is a closed curve, indeed an ellipse. In a curved space with Q^O, one 
has ^TZ, and consequently u^Q. The slightly non-circular orbit would not be 

a closed curve. It could be represented by a precessing ellipse, with two consecutive 
perihelia shifted by 

\ 1/2] 



dn r 
'dfTV 



(34) 



A0 = r(l] - w) = 27r(l - ^) = 27r 1 
where T = 27r/r2 is the orbital period. 
6. A Newtonian experiment 

Newtonian dynamics allows one to measure the circumferential f and curvature TZ 
radii of circular orbits by measuring the gravitational ac and centrifugal ac radial 
accelerations for a circular orbit, 

GM V' 

ac = ac = — . (35) 

7c 

In the Schwarzschild metric, the acceleration of a particle (a "planet" ) moving with the 
orbital velocity v along a circular orbit equals, 

a. = V.vl/ + V'^. (36) 

Here V = v/(l — t;^)^/^, and the scalars \1/ and R are expressed in terms of the time- 
symmetry Killing vector 77*, and the axial-symmetry Killing vector 

^ = -lHv'v^), B' = -%^. (37) 
In Schwarzschild coordinates this is, at the "equatorial plane" 6 = 7i/2, 

(rj^r^,)=gu = l-—, (^^6) = -r^ (38) 
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This allows one to calculate the results of the Newtonian experiment to measure the 
gravitational and centrifugal accelerations, 



In 



2M/r 



By comparing fl5S]) and fl5^ . one concludes that, 



r r 



r (1 - 2M/r 



,1/2 



7^(r) 



2M/r 



(39) 



(40) 



1 - 3M/r' 

Note, that to 0^{M/r) accuracy this is r = f = 7^. Therefore, curvature effects may 
appear at this order. We can also usefully calculate the derivative dr^{r)/dr, as only the 
derivative, not the absolute value of r*(r) is of interest. The following equation follows 
from the definition of Frenet's curvature radius TZ 

dr^ IZdf r — M 

dr r dr r — 3M 
The above formula allows one to write the metric of the 2-D space geometry, ds"^ 
dri + r'^d(j)'^, measured in this Newtonian experiment, 

2 



(41) 



ds' 



r - M 



dr"^ + rM 1 



r — 3M J \ r 

Inserting fHOj) into the Newtonian perihelion advance formula fl34|) one gets. 



-6Mr3 + 34M2r2 - Q2M^r + 3QM^ 



1 



(42) 



(43) 



Expanding this to the desired accuracy 0'^{M/r), one finally gets the same value for 
the perihelion advance as calculated in Einstein's theory. 



A0 



6n— + CM — 

r \ r 



(44) 



7. Light bending 

Knowing the space geometry, given by equation fH2|) . we may calculate the effect of light 
bending assuming that light travels along geodesic lines in space. In Newton's theory 
this assumption is equivalent to the Fermat principle, i.e. that light travels (with a 
constant speed) between two points A, B in space, minimizing the time travel Tab- 
The equation of motion for the coordinate is, in these circumstances, 
d^ r^A^drd^^ 

ds^ r(r - 2M) ds ds ' ^ ' 

from which we find ^ to be equal to 
d(h C 



ds r{r-2M)' ^^^^ 
The integration constant can be evaluated at the perihelion location r = Rq (i.e. where 
d(l)/ds = 0), yielding 

d(j) _ ^/RoiRo - 2M) 
'ds ~ r(r - 2M) 



(47) 



Using equations and (H7|) we find also 



dr _ r-m RojRo - 2M) 

ds r- M \l r(r - 2M) ' ^ 

After dividing equation ( 147|) by equation (148|) and substituting x = Ro/r the d(j)/dr 
equation can be integrated from i?o to oo (or x from to 1), wliicli will give us the half 
of n + 6. Let us also define fi = M/Rq, then 



;i - 2x/i)(l - 3x/i) y 1 - x^{l - 2/i)/(l - 2x/i) 
This integration can be expanded in a Taylor series for /i: 

77 + 6 dx (3x2 _^ 33. _ ^ 



yl — x^ Jo y/l — x^lx + l) 
As the first component on the right hand side is equal to 7r/2, we conclude that 

Jo V 1 - a; (a; + 1) -n^o 
Similar calculations in the Schwarzschild spacetime geometry give the same resuhl§| 

^«2./^^-^^4-. (52) 

Once again, the prediction of the Newtonian theory in the non-fiat space is found to be 
consistent with observations (and with Einstein's general relativity). 



8. Enhanced Newtonian Gravitational Theory 

Jiirgen Ehlers pointed out in 1961 that in Einstein's theory the curvature of the rest- 
space of irrotational matter is determined by its distribution and relative motion (see 
1221 in his article [3]). The equations governing such 3-space curvature for arbitrary 
irrotational flows are given in see their equation (54). Consequently it makes sense 
to consider gravitational dynamics in the context of 3-dimensional curved Riemannian 
spaces. As Newtonian theory is an approximation to General Relativity Theory, it is 
therefore interesting to see what happens in the case of Newtonian theory in a curved 
3-dimensional background space. 

In the case of isometric flows, 6 = = and there is a potential such that ii'^ = U^a 
where the gravitational potential U relates the Killing vector C, to the unit 4-velocity u"' 
by C,"" = e^u°' (see 1234 in Ehlers [3]). Then the relevant equation becomes 

^Rab = V, Vfef/ + VaUVbU +'^kp hab, (53) 

where Va is the 3-dimensional covariant derivative, p is the energy density of matter, 
and we have assumed anisotropic stress is zero {nab = 0) and a vanishing cosmological 



§ Which is twice the well-known flat-space and massive photon Newtonian prediction. 
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constant. Here hab = gab — UaUb is the metric of the three-spaces orthogonal to m". 
This case will include static spherically symmetric spacetimes. Taking the trace of this 
equation gives (see equation (55) in jl]) 

^R = 2kp, (54) 

where the potential terms have gone because of the relation between the 3-dimensional 
and 4-dimensional covariant derivatives. Together with the Poisson equation and 
equation of motion it defines the Enhanced Newtonian Gravitational Theory, 

'g'" V.Vfc<l> = -4nGp, 
Fi = mai, 

^R = 2kp. (55) 
For spherically symmetric spaces, the most general metric has the form, 

ds^ = A{r) dr^ + {dO^ + sin^ edcf)"^) , (56) 
and in the vacuum case, p = = one has 
3 ^' 

^rr = (57) 

'Roo = - I + l) si^' ^' (58) 

rA' 1 

'R,, = ^ - ^ + 1' (59) 
^ 2 2 2A' 

^ =^-I^ + Iv = °- 

Here the prime denotes a derivative with respect to r. Equation (!60|l has a unique 
solution, 

A(r)= (61) 

with tq being an integration constant. Its value cannot be determined by equations 
(135]) . but instead must be chosen by correspondence with experimentliJ]. Using the same 
procedure as in Sections [6] and U\ one proves that the choice tq = 4M gives the correct 
values for the perihelion advance and light bending (with accuracy (9(ro/r)). 

9. The two metrics 

We have shown that "experimentally" established and the "theoretically" postulated 
Newtonian metrics of the curved 3-D space corresponding to a spherically symmetric 
body are, respectively, 

11 In Einstein's theory, when one derives the Schwarzschild metric, a constant of integration is 
determined in a similar way, i.e by correspondence with Newton's theory. 
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(62) 

ds^ = ^1 _ dr' + [dO^ + sin^ Odcj)'') . (63) 

We have also shown that any spherically symmetric metric that obeys ^i? = must be 
isometric with f l63|) . The Ricci scalar for the "experimental" metric may be calculated 
to be 

\ ^)'^' r2[4(M/r)3 - 8(M/r)2 + 5(M/r) - 1] 

= 10{M/r)^ + 32{M/rf + ... 

= + O\M/r). (64) 
On the other hand a metric, 

has the Ricci tensor, 

CR)M^ = 2a(—] +4/3 (—] +...=0 + O\M/r). (66) 



Thus, the experimental metric (162!) and the theoretical metric (163|) describe, with 
accuracy (9^(M/r), i/ie same geometry of space. 



10. Conclusions 



We demonstrated that a Newtonian physicist may experimentally determine the 
geometry of the 3-D space ^gj^ by measuring gravitational and centrifugal accelerations. 
He may then predict by calculations the perihelion advance and the light bending as 
effects of the curvature of space. The predicted values agree with the ones measured. 
We also demonstrated that one may extend Newton's theory of gravitation by adding an 
equation that links Ricci curvature of space with the density of matter. We calculated 
the resulting theoretical metric of space ^gfj^ assuming spherical symmetry. In this 
metric, the values of perihelion advance and light bending also agree with those observed. 
The two metrics represent the same geometry, ^gj^ = ^gj^^ with accuracy (9^(M/r). 

Abramowicz [1] has shown that for spaces with constant Gaussian curvature Newton's 
theory predicts no perihelion advance. We speculate that this is why Gauss (and 
other XIX century mathematicians) who might have calculated Newtonian orbits in 
curved spaces, would have missed the effect of perihelion advance. Most probably, they 
would calculate orbits in spaces with a constant Gaussian curvature first. Gauss almost 
certainly made this calculation. He was a master in calculating orbits. He made himself 
famous at the age of 23 by calculating the orbit of Ceres, discovered in 1801 by Piazzi. 
He seriously considered the possibility that our space is curved. He even attempted to 
determine the curvature of space by measuring angles in a big triangle (69 km, 84 km, 
106 km) made by the summits of Brocken, Hoher Hagen and GroBer Inselsberg. Gauss 
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was not quick in publishing his results concerning curved spaces. It is known that he 
discovered most of Bolyai's results, but never published them. Gauss died in February 
1885, four years before Le Verrier discovered the effect of the perihelion of Mercury 
advance. 
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